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(Abstract  ) 

The  method  or  orthoqon alizeo  plane  waves  is 
applied  to  a  calculation  or  electronic  energy  levels 

IN  POTASSIUM,  USING  a  POTENTIAL  OBTAINED  r ROM  A  SELF- 

oonsistent  eielo.  The  energies  or  twenty-pour  states  at 

POUR  SYMMETRY  POINTS  IN  THE  BRILLOUIN  ZONE  HAVE  BEEN  OB¬ 
TAINED*  The  lowest  bano  is  dismissed  in  detail  ano  the 
qualitative  PEATUR C8  or  THE  density  or  STATES  are  pre¬ 
sented*  Departures  from  free  electron  banos  are  found* 
Comparison  or  higher  bands  with  slue  other  calculations 
suggests  that  certain  features  of  bano  schemes  may  be 
reasonably  independent  or  the  potential  useo* 


^Supported  by  the  Office  or  Naval  Research 


/ 


1.  INTRODUCTION 


Electronic  energy  levels  in  the  alkali  metals  h/'c 

BEEN  STUD  1 E  0  FOR  MANY  YEARS*  THE  MAJORITY  OF  THE  WORK 
HAS  CONCTRNEO  THE  LIGHTER  ELEMENTS,  LITHIUM  AND  SODILM^^ 


1.  An  extensive  bibliography  is  given  by 

J*  S*  Slater,  Technical  Report  No.  4  or  thi 
Solid  State  ano  Molecular  Theory  Group, 
f/.l.T.,  0953)*  UNPUBLISHED. 

Principal  attention  has  been  devoted  to  the  cohibive 

ENERGY  ANO  TO  RELATEO  PROPERTIES  SUCH  AS  THE  EQUILIBRIUM 
LATTIOE  CONSTANY  )t  *0  THE  COMPRESSIBILITY.  THE  HEAVIER  AL¬ 
KALI  METALS!  POTASSIUM,  RUBIOIUM,  ANO  CESIUM  HAVE  QEfN 

(2  -  7) 

8TUDIE0  LESS  EXTENSIVELY'  1 


2.  E.  Gorin,  Phys.  Zeits.  Sovjetunion  328 

3.  T.  S.  Ku  in  and  J.  H.  Van  Vleck,  Phys.  Rev,  72., 

3^2  0950). 

H.  Brooks,  Phys.  Rev.  21#  1027  (1955)* 

5-  Berman,  Callaway,  and  Wcoos,  Phys.  Rev.  101 , 

D467  0956). 

6.  F.  S.  Ham,  Solio  State  Physics  1^,  127  (1955)* 

7.  R.  M.  Sternheimer,  Phys.  Rev.  78*  235  (195^). 


THIS  WORK  H  C  RE  REPORTED  CONCERNS  AN  APPLICATION  O1'  TJ-E 

(Q) 

METHOO  OP  ORTHOGONAL I  ZED  PLANE  WAVES  TO  A  CALCULATION 


8.  C.  Herring,  Phys.  Rev,  £2,  1 1 69*  (1940)® 


Or  ELECTRONIC  ENERGY  BANDS  IN  POT A3S I UM • 

Gorin  made  a  calculation  of  the  cohesive  energy  of 

(?.}  (9  -  10)  AND  A 

POTASSIUM  USING  THE  CELLULAR  METHOO  ' 

9®  Eo  WlGNER  AND  F.  SEITZ,  PHYS,  REV., 

804  (1933)® 

10,  F«  C.  Von  oer  Lage  and  H»  A,  Bethe,  Phys,  Rev, 

21,  612  (1947)® 


POTENTIAL  OBTAINED  FROM  A  SELF-CONSISTENT  FIELD  FOR  THE 

•f 

k  ion.  The  results  of  his  work  were  in  markeo  disagreement 

WITH  EXPERIMENT,  IN  THAT  TOO  LITTLE  BINDING  WAS  OBTAINED, 

This  failure  gave  rise  to  the  belief  that  an  energy  ieve>. 
calculation  based  on  a  self-consistent  fielo  would  be  ve»y 
inacourate  for  the  heavier  alkali  met. lb.  The  quantim  or- 
FECT  METHOO  WAS  DEVELOPED  BY  KUHN  AND  VAN  VLEOKw/  AND  ST 
rlROOKS^  TO  AVOIO  THIS  DIFFICULTY  BY  USING  OBSERVEO  SPECTRO¬ 
SCOPIC  OATA  TO  CIROUMVENT  THE  OON8TRUCT  ION  OF  AN  EXPLICIT 

POTENTIAL,  IT  HA8  RECENTLY  BEEN  SHOWN,  HOWEVER,  9Y  BERKUN, 

(9  ) 

Callaway  and  Woods  '  that  if  proper  account  is  taken  of 
EXCHANGE  INTERACTIONS,  A  SELF-CONSISTENT  FIELD  METHOD  OAN 


GIVE  A  RESULT  OF  REASONABLE  ACCURACY  FOR  THE  COHESIVE  ENERGY 


OF  PC  TANS  #'JU*  n  Th  EM  BECOMES  INTERESTING  TO  EXTEND  THCIT 
WORK  TO  A  CALCULATION  OF  HIGHER  E LEOTRON I C  STATES* 

i 

In  making  this  EXTENSION,  the  method  of  orthcc.cm.vli2ED 
PLANE  WAVES  (OPW)  WAS  CHOSEN  IN  PREFERENCE  TO  THE  CELLULAR 
METHOO  BECAUSE  IN  THE  OPW  METHOD  IT  IS  NOT  NECESSARY  TO 
8ATISFY  BOUNDARY  CONDITIONS  EXPLICITLY*  THE  »ROBLSN  or 
ACCURATELY  SATISFYING  THE  BOUNDARY  CONDITIONS  IS  QU|VE 
S-tR  1 0U8  ANO  JlFFICULT  IN  THE  CELLULAR  METHOD'  ,  AMI  IT 

11,  F.  S.  Ham,  Ph.D*  Thesis,  Harvard  University 
(195^-)  (unpublished). 

WAS  FELT  THAT  WITH  THE  CELLULAR  METHOD,  IT  WOULD  HAVE  BEEN 

IMPRACTICAL  TO  STUOY  MORE  THAN  A  SMALL  NUMBER  OF  ENERGY 

STATES,  The  OPW  methoo  Has  THE  DISADVANTAGE  THAT  AN  EXPLICIT 

POTENTIAL  MUST  BE  USED,  AMO  MORE  IMPORTANT,  THAT  THE  ELECTRON 

(12) 

STATES  IN  THE  ATOMIC  CORE  MUST  BE  KNOWN'  .  I  HE  RE  A, RE 

12,  J,  Callaway,  Phys.  Rev.  22*  953*  (1955)« 

REASONS  WHY  THE  CtP.W*  METHOD  SHOULD  BE  EXPEOTEO  TO  WORK 
WEfcL  FOR  METALLIC  POTASSIUM*  AS  IN  THE  OTHER  AL'CAV  ME*AL'I. 
ONE  HAS  ONE  ELECTRON  MOVING  OUTStOE  OF  A  RELATIVELY  COMPACT 
ion  core.  Consequently,  overlaping  of  oore  wave  functions 
WILL  BE  NEGLIGIBLE*  THE  f.  ELF-CONS  I  STENCT  PROBLEM,  ALSO, 
SHOULD  NOT  BE  8ERIOUS  SINCE  IT  IS  REASONABLE  TO  EXPECT  THAT 
THE  AVERAGE  DISTRIBUTION  OF  ELECTRONS  IN  THE  CORE  SHOUL0  NOT 
BE  SIGNIFICANTLY  AFFECTED  BY  THE  VALENCE  ELECTRON*  THIS  1$ 


8UPP0  n  t.'O  8  V  A  CALCULATION  OF  MaRTREC  ANO  HARTREE  i-GH  THE  K 

IJOM  t  :i  WHICH  IT  WAS  FOUND  ''HAT  THE  WA Vi  FUNCTION  Or  ',lit  J 

+  f  \  X  i 

ELECTRON  DIFFERED  ONLY  SLVGHTLV  FROM  THAT  FOR  THE  !<  ON  '  '# 


M*  Oj  R.  HART  PEE  AN  '3  W.  HARTREE,  PROC,  CAM  3  T  t  !'  !  I. 

Phil.  Soc.  *£.,  550  (1936). 

» 

The  approximation  or  Wigner  and  Seitz  is  also  made;  that 

THERE  IS  ONLY  CNE  VALENCE  ELECTRON  IN  AN  ATOMIC  Of 'Li  UH  I CH 
MOVES  IN  THE  POTENTIAL  OF  THE  POSITIVE  ION,  ALL  OTHi;  CELLS 
BE  I  NO  NEUTRAL.  THE  CRYSTAL  POTENTIAL  AND  CORE  WAVE  TUMCTJCMS 
USEO  HERE  WERE  OBTAINED  FROM  A  SELF-CONSISTENT  FIELD  WITH 
EXCHANGE  FOR  THE  K  CON-  ^  . 


Vi,,  D.  R.  Hartrfe  ano  W,  Hartrec,  Proc,  Royal 
Society,  1o6a,  I4.50  (1956). 

Because  of  the  simplicity  of  the  physical  situm  on 

POTASSIUM  SHOULD  BE  A  VERY  FAVORABLE  SYSTELl'  FOR  AN  I  «■.:««» 

BANO  STUDY.  SUCH  A  STUOY  SHOULD  FURNISH  INFORMATION  CON¬ 
CERNING  •;  (1)  Validitt  of  thf  FREE  ELECTRON  APPROXIM, 

"i  H  E  R  F  ARE  THREE  ASPECTS  HERE  S  (a)  DEPARTURES  OF  EFFl  '  VI 
MAS  3  VALUES  FROM  (b)  PRESENCE  OF  ANO  HIGHER  T»  13  IN 

THE  E(k)  EXPANSION  ANO  (c)  DEPARTURES  OF  ENERGY  3URI  .  :  :fl 
FROM  SPHERICAL  SYMMETRY.  ONE  MIGHT  EXPECT  THAT  THM  .  EIMIC  S 
WO;- 4.0  r  :  MORE  PRONOUNCED  FOR  EXCITED  LEVELS  RATHER  T*  ui  FOR 
THE  LOWEST  BAND.  (2)  COMPARISON  OF  ENERGY  BANDS  IN  POTASSIUM 
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WITH  THOSE  OA..CULATED  FOR  OTHER  ELEMENTS  HAVING  THE  BOLT 
CENTERED  OU8IC  STRUCTURE?  ONE  OF  THE  MOST  IMPORTANT  PRO¬ 
BLEMS  IN  ENERST  BAND  THEORY  IS  THE  QUESTION  OF  HOW  SENSI¬ 
TIVE  IS  THE  BAND  STRUCTURE  TO  DETAILS  OF  THE  POTENT  s  AL* 

Comparison  of  energy  bands  in  different  elements  of  the  same 

CRYSTAL  STRUCTURE  SHOULD  FURNISH  INFORMATION  ON  THIS  QUESTION. 

In  aouition,  the  effect  of  ohanges  in  the  potential  on  the 

SANDS  IN  POTASSIUM  CAN  ALSO  BE  EXAMINEO. 

Unfortunately,  there  is  a  oearth  of  accurate  experimental 

EVIDENCE  PERTAINING  TO  ENERGY  BANDS  IN  POTASSIUM,  So  THAT 
A  DETAILED  COMPARISON  OF  THEORY  AND  EXPERIMENT  IS  NOT  POSSIBLE. 

Qualitative  comparisons  ano  pi;  eductions  can  be  made  in  some 
cases.  It  is  hopeo  this  theoretical  discussion  will  stimulate 
EXPERIMENTAL  EFFORT. 

A 

II.  THE.  CRYSTAL  POTENTIAL 

We  ASSUME  TiAT  EACH  ELECTRON  EXPERIENCES  ONLY  THE 
POTENTIAL  OF  THE  C OR  £ S POND  I  NO  POSITIVE  ION.  THIS  POTEN¬ 
TIAL  CONSISTS  OF  TWO  PARTS?  (1)  THE  COULOMB  POTENTIAL  OF 
THE  ATOMIC  NUCLEUS  ANO  THE  AVERAGE  DISTRIBUTION  OF  THE  CORE 
ELECTRONS,  AND  (2)  THE  EXCHANGE  INTERACTION  BETWEEN  THE  V AL*" 
ENCC  ELECTRONS  ANO  THE  CORE.  THE  COULOMB  POTENTIAL  IS  OBTAINED 
IN  THE  OBVIOUS  WAY  FROM  THE  CORE  ELECTRON  0 1 8TR I  BUT  ION.  THI 
EXCHANGE  INTERACTION  18  APPARENTLY  QUITE  IMPORTANT  IN  OBTAINING 
NUMERICAL  RESULTS  IN  A  COHESIVE  ENERGY  CALCULATION. 
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We  CAN  DEFINE  AN  EFFECTIVE  EXCHANGE  POTENTIAL  F OP  THE 
•  TATE  yj  (  I  *  TH  IRRCOIIC  ISLE  REPRESENTATION  OF  WAVE  VECTOR  K_) 
IN  THE  F 0  LLOW I N  G  WAY? 


(  «  ,  )* 


NLM  ^  *2  ^ 


Ri2  a 


(*2) 


*D 


W  (r  ) 

tnlm'  i' 


* 

n 


H 


WHERE  t  REPRESENTS  A  CORE  SAVE  F  UNOT  (t  ON «  THE  EXCHANGE  POT«'»‘IIAL 
H  LM 

COMPUTED  ACCORDING  TO  (1)  WILL  OIFFER  FROM  STATE  TO  STATE, 

AND  OBVIOUSLY  DEPENDS  ON  THE  WAVE  FUNCTION  OF  THE  STATE  CON* 
SIOEREO.  Th 1 8  MEANS  THAT  A  SELF-CONSISTENT  SOLUTION  Of  THE 
HARTREE-FOOK  EQUATION  IS  REQUIRED.  SLATER  HAS  PROPOSCO  TWO 

(15) 

METHOOS  OF  AVERAGING  THE  EXCHANGE  POTENTIAL  " 

15«  J*  C.  Slater,  Phys.  Rev.,  til^,  5^5  (1951)* 


WHIOH  ALLOW  ONE  TO  USE  ONE  EXCHANGE  POTENTIAL  FOR  ALL  STATES. 


IT  18  OOUBTFUL  THAT  8U0H  A  PROCEDURE  WILL  HAVE  QUANTITATIVE 


3UOOESS 


(16),  (17) 


16.  Herman,  Callaway,  and  Aoton,  Phys.  Rev.  95. 

571,  (5954). 

17*  J.  Callaway,  Phys.  Rev.,  500  (1955)* 

A  more  accurate  proceoure  woulo  be  to  assume  the  exchange 

POTENTIAL  FOR  A  GIVEN  STATE  OEPENQS  PRIMARILY  ON  THE  ANGULAR 
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MOMENTUM  OF  THE  STATE  CONSIDERED,  OR  IN  THE  CASE  OF  THE  SOLID, 

ON  THE  PREDOMINANT  ANGULAR  MOMENTUM  IN  THE  OE  0CI..P08 1 T  ION  Or 

S*  INTO  SPHERICAL  HARMON  I CS ^ ^ ^ .  |n  OBTAINING  THE  EXCHANGE 
POTENTIAL  TOR  A  STATE  OF  ANGULAR  MOMENTUM  L,  OF  A  VALENOE 
ELECTRON,  IT  IS  QUESTIONABLE  WHETHER  IT  IS  DESIRABLE  TO 
AVERAGE  THE  EXCHANSE  POTENTIALS  OF  THE  CORE  ELECTRON  STATES 
OF  ANGULAR  MOMENTUM  L,  8INCE  THESE  HAVE  ENERGY  VALUES  IN  GENERAL 
FAR  BELOW  THAT  OF  THE  STATE  WE  ARE  CCNSIDERIN3.  IT  SEEMS  MORE 
REASONABLE  TO  CONSTRUCT  AN  EXCHANGE  POTENTIAL  FOR  A  GIVEN  L 
FROM  AN  APPROXIMATE  WAVE  FUNCTION  FOR  A  VALENCE  ELECTRON 
STATE  OF  THAT  L. 

Such  a  procedure  was  followed  in  this  caloulat it n*  An 
s  state  exchange  potential  was  constructed  from  the  self* 

C<  NS  1  STENT  riCLO  CORE  wave  functions  and  the  lowest  ortho- 
GONALIZEO  PLANE  WAVE  FCR  THE  STATE  f  .  THIS  POTENTIAL  WAS  USEO 
FOR  THE  STATES  f(,  H(,  Pjf  AND  N(.  IT  IS  TO  BE  NOTED  THAT 
THE  STATE  N  WILL  CONTAIN  AH  AOMIXTURE  OF  D  FUNCTIONS.  A  P 
STATE  EXCHANGE  POTENTIAL  WAS  SIMILARLY  OBTAINED  FROM  THE  LOWEST 
C.P.W.  FOR  THE  STATE  H^,  ANO  USEO  IN  THE  CALCULATIONS  FOR 
r^,  H^,  P^,  N('  Nj*,  N^*.  P|^  WILL  AL80  HAVE  SOME  D  CHARACTER. 

A  O  STATE  EXCHANGE  POTENTIAL  WAS  OBTAINEO  FROM  A  0  FUNCTION 
OALCULATEO  WITHOUT  EXCHANGE  IN  A  PRIOR  CALCULATION  OF  TERMS 

IN  THE  E  (±)  RELATION*  THIS  POTENTIAL  WAS  USED  FOR  THE  O-LIKE 
STATES  r^y  ri2,  H Zy  H12,  Py  Ng,  N y  Nj^,  AND  AL 80  FOR  THE  F 
LIKE  STATES  r^y  f*p  8  ,  H ^y  Hg  I ,  P^  ANO  Ng  I  *  IT  WAS  ASSUMED  IN 
USING  THE  O.P.W.  METHOD  THAT  THE  CORE  STATES  03UL0  WITH 


SUrriOtENT  ACCURACY ,  BE  CONSIDERED  A8  EIGENFUNCTION 
IN  THESE  POTENTIALS. 

One  unsatisfactory  feature  of  this  procedure  is 

that  in  THE  C A8E  OF  THE  S  AND  P  STATE8  THE  EXCHANQE 

POTENTIAL  AOOORDINS  TO  (1)  WILL  HAVE  INFINITIES  VfHERE 

THE  APPROXIMATE  <| T‘  HAS  ZEROS.  THESE  WERE  REMOVED  IN  SUCH 

K 

A  WAY  AS  70  OIVTC  A  REASONABLY  SMOOTH  EXCHANQE  POTENTIAL* 

The  coulomb  potential  At.o  the  potential  for  s,  p,  ano 
o  states  are  tabulated  in  Table  t.  Certain  irregularities 

MAY  BE  NOTED  IN  THE  QUANTITIES  RV  AND  RV  IN  THE  NEIGH* 

S  P 

B0RH00D  OF  R1^.  THESE  RESULT  FROM  THE  METHOD  OF  TREATING 
THE  EXCHANQE  POTENTIAL  NEAR  A  SINGULARITY.  SINCE  THESE 
IRREGULARITIES  OOOUR  IN  A  REGION  IN  WHICH  THE  WAVE  FUNCTION 
18  SMALL,  THE  EFFECT  ON  THE  ENERGY  SHOULD  NOT  BE  LARGE. 


III.  CALCULAT ION  OF  THl  ENfcRGY  LEVELS 
The  O.P.W.  mcthoo  leads  to  a  secular  equation  or 


THE  FORM 


OET  (XK,  HXh)  -  E(X„,  XM)  -  0 


WHCRE44  18  THE  ORTSTAL  HAMILTONIAN  ANO  XK  IS  AN 
ORTHOGONAL IZEOv PLANE  WAVE! 

xk  ■  “  'hi  “«  / -'£v  *»  (5) 

Here  is  a  oore  function  for  the  core  state  u,  Rg 

IS  A  LATTIOE  VE0T0R  ANO  *»„■  ■■  ■■  /3  /4 .  (ft)  C  dff. 

WITHV^q  *  THE  VOLUME  OF  THE  ATOMIC  CELL  (Ivi). 
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Ido  This  calculation  was  made  for  a  lattice 

PARAMETER  OF  5«20A  ,  OR  ATOMIC 

UNITS# 


In  PRACTICE  WE  USE  LINEAR  COMBINATIONS  C(  0*tHG QOMAL IZ ED 
PLANE  WAVES  WHICH  TRANSFORM  ACCORDING  TO  PARTICULAR 
IRREDUCIBLE  REPRESENTATIONS  IN  SETTING  UP  (2). 

We  assume  that  the  core  states  4.  arc  eigenfunctions 

J 

OF  THE  CRYSTAL  HAMILTONIAN#  IN  THIS  CASE,  THE  MATRIX 
ELEMENTS  in  (2)  are 

(VHXh,',Z£*  *v  ^  -  ZJCJ  “'kj'S.j  Vi) 


ANO 


^XK*  XH^  "  *Scj 


WHERE  £  IS  THE  ENERGY  VALUE  OF  THE  CORE  STATE  J  ANO 

V(k)  is  a  Fourier  coefficient  of  potential 


V(K)  A 


I  K  •  R  V 


<«) 


(5) 


The  Fourier  coefficients  V^k j  are  siven  in  a  table  (2) 

2  a  2 

AS  FUNCTIONS  OF  M  ■(  -g-  «)  FOR  S,  P  ANO  0  STATES  (a  18  THE 

lattice  parameter)# 

In  connection  with  the  Fourier  coefficients.  Butcher  (1J) 


15c  Pc  No  PliTOHER,  PROC.  RHYS,  SOO#  6Ua.  763  095U 


HAS  OBTAINED  V(m*“2)  FOR  POTASSIUM  ANO  SEVERAL  OF  THE 


OTHER  ALKALI  METALS  FROM  A  DISCUSSION  OF  THE  OPTICAL  PROPERTIEB 
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8ASE0  ON  THE  APPROXIMATION  OF  NEARLY  FREE  ELECTRONS •  H  18 
VALUE  OF  0*305  E.V.  SEEMS  TO  BE  TOO  SVALL  BY  NEARLY  AN 
ORDER  OF  MAGNITUDE.  THIS  IS  PROBABLY  DUE  TO  FAILURE  OF  THE 
APPROXIMATION  OF  NEARLY  FREE  ELECTRONS. 

r 

O.P.W.  EXPANSIONS  WERE  CONSTRUCTED  FOR  2i|.  STATE8  MEN* 
tioneo  in  Sect.  II  at  the  symmetry  points  f“,  H,  P,  and  N  in 
the  Brillouin  Zone.  The  Brillouin  Zone  for  the  body  centered 

CUBIC  LATTICE  IS  3H0WN  IN  FIGURE  I.  LINEAR  COMBINATIONS  OF 
ORTHOGONAL! ZED  PLANE  WAVES  WERE  EMPLOYED  WHICH  TRANSFORMED 
A0C0R0IN6  TO  THE  IRREOUCISLE  REPRESENTATIONS  OF  INTEREST 
(SEE  APPENOJX  OF  REFERENCE  17).  BECAUSE  HIGH  SPEED  ELECTRONIC 
COMPUTING  EQUIPMENT  IS  NOT  AVAILABLE  AT  THE  UNIVERSITY  OF 

Miami,  the  computations  010  not  in  general  involve  higher  than 

FOURTH  ORDER  DETERMINANTS.  FIFTH  CROER  DETERMINANTS  WERE 

Solve 0  for  twc  states:  A«o  P^«  Third  and  secono  order 

DETERMINANTS  WERE  USED  FOR  THE  STATES  IN  THE  F  BANO  WHERE  ONLY 
A  QUALITATIVE  INDICATION  OF  THE  ENERGY  IS  REQUIRED.  SUOH 
HIGHLY  EXCITEO  STATES  CAN  PROBABLY  BE  REASONABLY  WELL  REPRE¬ 
SENTED  BY  SMALL  NUMBERS  OF  PLANE  WAVES.  THE  LOWEST  TWO 
EIGENVALUES,  (IN  ONE  CASE,  THREE)  OF  ALL  THE  REPRESENTATIONS 

A 

CALCULATED  ARE  GIVEN  IN  TABLE  3«  The  ORDER  OF  THE  LEVEL8 
AT  THE  FOUR  SYMMETRY  POINTS  IS  SHOWN  IN  FIGURE  2. 

In  croer  to  stuov  the  convergence  of  the  O.'.W.  ex¬ 
pansion,  AN  eighth  order  seoular  determinant  for  the  STATE  r 

•A3  8OLVE0  ON  THE  WHIRLWIND  COMPUTER  AT  THE  MASSACHUSETTS 
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Institute:  of  Technology  through  the  courtesy  of  Professor 
J.  Cc  Slater  and  Mr.  F.  J.  Corbato.  Convergence  of  the 

EIGENVALUE  TO  THE  VALLE  “.iU|-7  FOUNO  BY  THE  CELLULAR  METHOD 

FOR  THIS.  STATE  FOR  THE  POTENTIAL  V  APPEARED  TO  BE  SLOW. 

S 

(The  BEST  VALUE  OS'  THE  ENERGY  OF  THE  LOWEST  STATE  IS  *olj-55 
RY0BERG8  FOUNO  FROM  THE  HARTREE-FOCK  EQUATIONS  DIRECTLY 
WITHOUT  THE  APPROXIMATION  OF  AN  EXCHANGE  POTENTIAL.)  TABLE 
SHOWS  THE  LOWEST  EIGENVALUE  IN  EACH  ORDER:  FROM  FIRST  TO 
EIGHTH,  AND  THE  LOWEST  SIX  EIGENVALUES  IN  EIGHTH  OROE R .  IT 
IS  INTERESTING  TO  NOTE  THAT  THE  MAJOR  CHANGE  OCCURS  BETWEEN 
FIRST  ANO  SECOND  ORDER  (WHERE  THERE  ARE  13  WAVES)*  THERE 
ALSO  SEEMS  TO  8E  A  CLUSTERING  OF  EIGENVALUES  IN  A  REGION 
NEAR  £.*=+1023.  A  SMALL  EIGENVALUE  CHANGE  IN  EACH  ORDER  IS 
NECESSARY,  BUT  NOT  SUFFICIENT  TO  ENSURE  CONVERGENCE., 

Using  this  example  as  a  guide,  we  would  expect  that 

TH08E  SOLUTIONS  OF  FOURTH  ORDER  DETERMINANTS  BASED  ON  FORTY 
OR  MORE  WAVES  SHOULD  BE  CONVERGENT  TO  ABOUT  0 .  Oij.  RYDBERG. 

This  estimate  may  be  ccnserv at ? y»i,  for  some  of  the  higher 

STATES  MAY  BE  BETTER  REPRESENTED  BY  A  FEW  OR T HO GO N At  I Z E D 
PLANE  WAVES  THAN  THIS  STATE.  ?N  PARTICULAR,  THE  CONVERGENCE 
OF  THE  0  BAND  STATES  MAY  DE  StMEWHAT  BETTER.  IT  IS  ALSO 
REASONABLE  TO  EXPECT  THAT  THE  DIFFERENCES  BETWEEN  ENERGY 
LEVELS  FOUNO  U81NG  APPROXIMATELY  THE  SALE  NUMBER  OF  WAVES 

is  stable.  For  this  reason,  it  is  likely  that  the  state 

WILL  CONTINUE  TO  LIE  BELOW  N  f  '  „  HOWEVER,  IT  IS  POSSIBLE 
THAT  IN  HIGHER  ORDER,  Pj^  WOULD  BE  FOUNO  TO  LIE  BELOW  H^f 
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8  I  NOE  H,|j  APPEARS  NEARLY  0  ON  PER  61  NT  I  N  FOURTH  ORDER  WITH 
5U  WAVES  ,'i/H  I  LE  P.  J$  FARTHER  FROM  CONVERGENT  IN  FIFTH  ORDER 
WITH  c\\.  STAVES.  (By  A  P  E  CU  L  I  I.tt  I  ▼'*  OF  THE  GROUP  THEORY,  EACH 
ORDER  OF  J'HE  SEOUL  \ T  DETERMINANT  FOR  P^  ONLY  INCLUDES  A 
SMAU.  NUMBER  OF  ADDITIONAL  V,|,VES„  )  MANY  OF  THE  RESULTS  OF 
THIS  liORK  OEPENO  ONLY  ON  RELATIVE  POSITIONS  OF  THE  LEVELS? 

AND  At  E  PROBABLY  R  EL  1  ABLE  EVEN  #1*  THE  ABSOLUTE  VALUE  OF 
THE  ENERGIES  ARE  SDMEWHAT  UNCERTAIN#  NEVERTHELESS,  THE  CON¬ 
VERGENCE  APPEARS  TO  BE  MUCH  POORER  THAN  THAT  OF  THE  AUGMENTED 

(20) 

PLANE  WAVE  METHOD' 


20.  D.  J.  Howarth,  Phvs.  Rev#  ^2.,  U69  (1955)< 


A  VERY  IMPORTANT  QUESTION  IN  TKE  CALCULATION  OF  ENERGY 
BANDS  IS  THE  SENSITIVITY  OF  THE  ELECTRON  ENERGY  LEVELS  TO 
CHAN  G 1  S  IN  THE  POTENTIAL*  ONE  WOUIC  HOPE  THAT  SUCH  FEATJRES 
AS  THE  RELATIVE  OROER  OF  THE  LEVELS  AT  A  SYMMETRY  POINT 
WOULD  BE  REASONABLY  INDEPENDENT  OF  DETAILS  OF  THE  POTENTIAL. 

Unfortunately,  the  contrary  result  seems  to  be  inoicateo  by 

(20  1 

SOME  WORK  OF  HoWARTh'  HOWARTH  FOUND  THAT  THE  ASSUMPTION 

THAT  THE  POTENTIAL  IS  CONSTANT  IN  THE  REGION  OUTSIDE  THE 

INSCOSBEO  SPHERE  IN  THE  ATOMIC  CELL  WAS  SUFFICIENT  TO  INVERT 

THE  TRIPLY  DEGENERATE  AND  DOUBLY  DEOENERATE  0  BAND  LEVELS  AY 

THE  OENVER  OF  THE  BRILLOUIN  ZONE  IN  COMPARISON  TO  AN  EARLIER 

(21  ) 

CALCULATION.'  THIS  MODIFICATION  OF  THE  POTENTIAL  IS  USED 


21  n  K.  J  *  Hcwarth,  Proc.  Roy.  Soc.  (Lonoon)  A22u 

513  0953  '  • 
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NOT  ONLY  IN  THE  AUGMENTED  PLANE  WAVE  METHOD  EMPLOYED  BY 
HOWARTH^1^  BUT  ALSO  IN  THE  METHOD  OP  KOHN  AND  ROSTOKER^*^* 

22.  W.  Kohn  and  N.  Rostoker,  Phys.  Rev.  9ll.  1111,  (195*+)» 

The  EFFECT  OF  SUCH  A  MODIFICATION  OF  THE  POTENTIAL  WAS  EXAMINED 
in  this  case.  Since  the  core  wave  functions  oo  not  extend 
INTO  THE  REGION  IN  WHICH  THE  POTENTIAL  IS  TO  DE  MCOIFIEO, 

IT  IS  NECESSARY  ONLY  TO  CONSIDER  THE  EFFECT  OF  APPROPRIATE 
OHANGES  IN  THE  FOURMER  COEFFICIENTS  OF  POTENTIAL.  To 

SUFFICIENT  ACCURACY,  WE  MAY  ASSUME  THAT  THE  ORIGINAL  POTENTIAL 

p 

18  JUST  THE  COULOMB  POTENTIAL  /  OF  THE  POSITIVE  ION.  IT  IS 

R 

ADVANTAGEOUS  TO  CHOOSE  THE  CONSTANT  POTENTIAL  So  THAT  THE 

P 

M  *0  COEFFICIENT  OF  THE  DIFFERENCE  VAN  1 8HES •  SINCE  THE  COULOMB 
POTENTIAL  IS  REASONABLY  FLAT  IN  THIS  REGION,  THE  0 1 FFEREN  E  IN 
THE  POTENTIALS  IS  SMALL,  ANO  THE  FOURIER  COEFFICIENTS  OF  THE 
DIFFERENCE  IN  POTENTIAL  (WHICH  MAY  BE  COMPUTED  ANALYTICALLY) 

ARE  SO  SMALL  THAT  THE  EFFECTS  ON  THE  ENERGY  LEVELS  ARE  QUITE 
NEGLIGIBLE.  THERE  MAY,  HOWEVER,  BE  GOOD  REASON  WHY  THIS 
MODIFICATION  IS  NOT  SO  SERIOUS  AS  IN  HOWARTH'S  CASE^°^* 

IV.  DEPARTURES  F  RO.V  THE  EFFECTIVE  il'.ASS  APPRl  X  lf,AT  ION 

The  effective  mass  near  k*o  can  se  taken  from  the  work 
of  Berman,  Callaway  ano  lvooos  as  per  oent  of  the  free 

CLEOTRON  MASS.  HOWEVER,  THE  PROXIMITT  OF  0  BAND  STATES  AT 
N,  H,  ANO  P  TO  THE  GROUND  STATE  T(  (NOTE  THAT  N(,  ANO  Pj, 

CONTAIN  0  FUNCTIONS  AS  WELL  AS  S  ANO  P  RESPECTIVELY)  SUGGEST 


THAT  THERE  MAY  3E  APPRECIABLE  K  TERMS  IN  THE  EXPANSION 

of  E  (kJ  near  £*0.  The  coefficient  of  oan  be  oalcu- 

LATEO  ACCOROING  TO  THE  PROCEDURE  OF  SlLVERMAN^2^  WITHIN 


23.  R.  A.  Silverman,  Phys.  Rev.  227  (1952). 


THE  FRAMEWORK  O'  THE  CELLULAR  METHOD. 

According  to  Silverman,  if  we  write 
E  (k)  «  E0+E2k2+E)+k^ 


THEN 


2  2 

V  V  %  £2  - 
5 


JlisS2  <<  A>" 

15  Y 


+  7  2  --- 

“7*?.  L  r  v>  s 


r  '  E,  f*  \  \  -  R.  W  («  )  /  P  ’2  ON 


3  O  S  O  P 


S'  0 


PDZ  (r.) 
P  s 


WHERE  R  IS  THE  RAOIUS  OF  THE  ATOMIC  SPHERE 
S  ^ 

f  *  -.V  <"s>  - 


U  is  THE  CELLULAR  METHOD  WAVE  FUNCTION  FOR  K-O,  /  IS 

o  0 

A  SOLUTION  OF  THE  WAVE  EQUATION  FOR  A  0  STATE  OF  ENERGY 

E  ,  ANO  P„  IS  THE  P  FUNCTION  INVOLVED  IN  THE  USUAL  CALCULA- 
O  P 

TION  OF  THE  EFFECTIVE  MASS.  THIS  EXPRESSION  WAS  EVALUATED 

V 

USING  THE  RESULTS  FOR  U„,  E  .  E„  ,  AND  P_  OF  THE  PREVIOUS 

O  O  *2  P 

CALCULATION  OF  BERMAN,  CALLAWAY  ANO  MfOOOS^^  PLUS  A 

0  FUNCTION  OALCULATEO  USING  THE  POTENTIAL  V.  GIVEN  IN  TABLE  I 

D 
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FOR  ENERGY  E^.  THE  QUANTITY  (  L °0 )  WAS  EVALUATEO 

0  W  V  eo 

BY  EXPRESSING  UQ  FOR  R>4  AS  A  SUM  OF  REGULAR  AND  IRREGULAR 
CO  ULOM3  FUNCTIONS  WITH  COEFFICIENTS  THAT  ARE  KNOWN  FUNCTIONS 
OF  ENERGY  ACCORDING  TO  THE  PROCEDURE  OF  HAM  AND  THEN 

24-  F.  S.  Ham,  ONR  Technical  Rfport  No.  204*  Ckuft 

Laboratc i?v,  H/.rvaro  University,  ( 1 95^4- )  (unpublisheo). 

DIFFERENTIATING  THE  EXPANSION.  The  VALUE  OF  E  *  P  T  A  1 IV  E  P  FAS 
1.64  *N  ATOM  1 0  UNITS.  SUCH  A  TERM  WOULD  GIVE  A  CONTRIBUTION 
TO  THE  COHESIVE  ENERGY  IN  THE  CELLULAR  METHOD  OF 

OR  5.4  K~CAL/tA0L  'A  REPULSION),  WHICH,  HOWEVER,  IS  PR08ABLT 
CANCELLED  BY  THE  ADDITIONAL  ATTRACTIONS  RESULTING  FROM 
POLARIZATION  AND  RELATIVISTIC  EFFECTS. 

This  relatively  large  value  of  Ej^  indicates  that  there 

WILL  BE  MARKEO  DEVIATIONS  OF  THE  ENERGY  SURFACES  FROM  THE 
FREE  ELECTRON  APPROXIMATION.  AS  POINTED  OUT  BY  DR.  M.  H. 

(2S  ) 

Cohen  ,  a  large  Ej^  also  suggests  that  there  may  be 
25«  M.  H.  Cohen,  Private  Communication. 

IMPORTANT  DEVIATIONS  O1-'  THE  ENERGY  SURFACES  FROM  SPHERICAL 
SYMMETRY,  SINCE  IF  EXP  AN  OED  IN  POWERS  OF  JC,  THE  LOWEST  8UCH 

I 

TERMS  WILL  BE  PROPORTIONAL  TO  THE  FOURTH  POWER  OF  THE 


WAVE  VECTOR 


V,  STRUCTUKt  OF  THE  LOWEST  f  AND 


In  order  to  determine  the  form  of  the  lowest  E  (j<) 

SURFACE,  IT  as  NECESSARY  TO  DETERMINE  THE  ENERGIES  OF 
A  LARGE  NUMBER  OF  THE  STATES  BETWEEN  THE  END  POINTS  f 

« 

AND  H„„  OF  THE  3ANO,  THIS  IS  / ERY  DIFFICULT  BECAUSE  A 

12 

VERY  LARGE  NUMBER  OF  TERMS  IN  THE  0„P,W.  SECULAR  DETERMINANT 
WOULD  BE  REQUIRED  IN  OROEF  TO  INCLUDE  A  REASONABLE  NUMBER 
OF  WAVES,  EVEN  ALONG  S  YMMf  TRY  AXES  LIKE  A(100).  ACCORDINGLY, 
WE  HAVE  RECOURSE  TO  AN  I  N  T  E  RPOl.  A  T  !  0  N  SCHEME.  IT  SEEM8 
NATURAL  TO  EXPAND  E ( K  )  IN  A  Fo  JR  I  ER  SERIES,  INCLUDING  ONLY 
THOSE  TERMS  WHICH  HAVE  TH  l  PR0»FR  SYMMETRY,,  AN  EXPANSION 
JN  POWERS  OF  K,  I  „  E  .  ,  KUEIO  H  AH  MO  J  II  CS  ABOUT  K»0  OOES  NOT 
SEEM  OES'RABLE  BECAUSE  IT  WILL  BE  DIFFICULT  TO  OBTAIN  THE 
PROPER  BEHAVIOR  Of  E(X’  AT  SYMMETRY  POINTS.  THE  FOURIER 
SERIES  APPROACH  MEETS  ’Hi'  REQUIREMENT.  If  WE  CONSIDER  THE 
LOWEST  PLANE  WAVES  THAT  ATE  PERIODIC  IN  THE  RECIPROCAL 
LATTICE,  WE  HA  VC 

E(jO»E0+  co  3j  oo*y  coif  «-3(oos2j‘  +coc  Z  J  *cos2j'  ) 

+Y(cos2  (  cos2/V->cos2  co$2f  -*"cos2  'fl  cos 2  £ 

/  '  (8) 

+  (cos  3^  'ZQIi"j  CO 

WHERE  j  *SKXA»  E  fc°  )• 

Such  an  expansion  will  ha  e  the  proper  zero  gradient  at 

SYMMETRY  POINTS,  THIS  IS,  IN  FACT,  JUST  THE  FORM  OF 
EXPRESSION  FOR  A  SINGLE  S  BAND  THAT  IS  OBTAINED  FROM  THE 
TIGHT  BINDING  I U T C R PC L A T  I  C N  SCHEME  OF  SLATER  AND  KOSTERV  ' 


rJi\  I  n  ,  ,,  r •  r  Incrro  P"WS  .  Rev.  Gi±_ 

cU  0  U  *  O  .J  ■  H  H  l.  r\  r»  •«  -*  *  .  v  nc*  *  0  V  .*y, » 


1o»j.8  ^9=4) 


An  attempt  was  made  to  I  |J  the  lowest  sand  with 

AN  EXPRESSION  0*  TM  S  'cRU.  A  f  E  VE  PARAMET-R  SCHEME 

r 

WAS  CHOSEN  USING  Eq  ,  ,  YT  d*>  AS  IN  (8).  THE  CON¬ 

STANTS  WERE  DETERMINED  FROM  THE  ENERGY  OF  THE  T  STAVE 
(Ep  »  AND  THE  EFFECTIVE  MASS  NEAR  i<  “0  TAKEN 

‘J 

FROM  THE  CALC ELAT  10  i  OF  8ERMA.N ,  CALLAWAY  AND  WOODS 

*1.l68)(5),  a;  o  the  lowest  energies  of  the 
sta'S^  H  „  P^t  amo  The  constants  are  given  in 

Table  V, 


IaB'J  V 

Parameters  in  tc  a  i  on  t  J  ; : 

E  -  -,2l:o 

r 

c  - 

r-  a  -*0277 

Y  *  ■*•  0106 

t  -  +.0i  2 

The  energy  bands  determimEO  by  equation  (3)  arc  shown 

IN  FIGURES  ‘1  AMD  l(.t  FOR  'HE  00  AND  110  AXES,  WHERE  THEY 
ARE  COMPARED  WITH  THE  PARABi  l  !C  BAUDS  FOR  M E REE  ELECTRONS 
WITH  AN  EFFECTIVE  MASS  RATIO  ‘  /fl*  e  1-168- 


I? 
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The  bands  shown  in  figures  3  a  no  4  deviate  in 

IMPORTANT  RESPECTS  FROM  THE  FREE  ELECTRON  APPROXIMATION. 

It  IS  LIKELY  THAT  THESE  DEVIATIONS  ARE  TOO  LARGE,  SINCE 
IF  £j^  IS  CALCULATED  FROM  (8),  IT  TURNS  OUT  TO  BE  CON¬ 
SIDERABLY  LARGER,  THOJGH  OF  THE  SAME  SIGN,  THAN  THAT  FOUND 
FROM  THE  CELLULAR  METHOO,  IT  IS  NOT  EASY  TO  FIT  A  SPHERICAL 
FREE  ELECTRON  3AN0  ACCURATELY  IN  THIS  WAY.  HOWEVER,  THE 
QUALITATIVE  FEATURES  MAY  BE  CORRECT.  THE  BANDS  FIRST 
RISE  SOMEWHAT  FASTER  THAN  THE  FREE  ELECTRON  APPROXIMATION 
INDICATES  ANO  THEN  LEVEL  OFF  SO  AS  TO  HAVE  ZERO  SLOPE 
COMING  INTO  H  ANO  N. 

IT  WAS  .  BEL  IE VEO  THAT  ACCURACY  OF  EQUATION  (8) 

WITH  THE  FIVE  PARAMETERS  :lVtN  WAS  NOT  SUFFICIENT  TO 
JUSTIFY  A  OETAILEO  AND  ACCURATE  CALCULATION  OF  THE  OENSITV 
OF  STATES.  THE  NUMERICAL  CALCULATION  OF  A  DENSITY  OF 
STATES  FROM  (<j)  WOULO  HAVE  TO  BE  DONE  VERY  CAREFULLY 
8INCE  WE  ARE  LOOK  I  M3  FOR  DEVIATION  FROM  THE 

BEHAVIOR  OF  A  FREE  ELECTRON  BAND.  IN  ORDER  TO  GET 
SOME  INSIGHT  INTO  THE  QUALITATIVE  FEATURES  OF  THE  OENSITY 
OF  STATES  CURVE,  IT  IS  PROBABLY  SUFFICIENTLY  ACCURATE  TO 
USE  THE  METHOD  OF  HOUSTON  (27)  IN  WHICH  THE  DEN  I  TV  OF 

27.  W.  Vo  Houston,  Rev,.  Woo.  Phys.  20,  l6l  (194^)» 


19 


STATE8  CUR VC  IS  COMPUTED  ALONG  CERTAIN  SYMMETRY  LINES 
IN  THE  BRILLOUIN  ZONE,  ANO  THE  OVER  ALL  DENSITY  OF  STATES 
8S  FOUND  BY  IN  TE  R  PO  L  A  T  C  N .  HOUSTON’S  METHOD  IS  KNOWN 
TO  GIVE  USE  TO  SPURIOUS  SINGULARITIES  RESULTING  FROM 
EXAGGER/TED  CONTRIBUTIONS  FROM  SYMMETRY  POINTS  SUCH 
AS  N,  H,  ANO  P  WHERE  THE  ENERGY  SURFACES  HAVE  ZERO 
3RADIENT.  BUT  THE  OCCUPIED  PORTION  OF  THE  BAND,  WITH 
WHICH  WE  ARE  PRIMARILY  CCNCT  RNED,  DOES  NOT  EXTEND  UP  TO 
N,  WHICH  IS  THE  CRITICAL  POINT  CLOSEST  TO  THE  ORIGIN. 

Consequently,  Houston’s  method  will  sive  no  spurious 

S  INGULAIII  T  |ES  IN  THE  OCCUPIED  FORT  ION  OF  THE  DENSITY  OF 
STATES  CURVE. 

A  QUALITATIVE  DENSITY  OF  STATES  CURVE  IS  SHOWN  IN 
FIGURE  5  WHERE  IT  IS  COMPARED  WITH  THAT  FOR  A  PARABOLfC 
BAND.  FOR  ENERGIES  ONLY  SLIGHTLY  ABOVE  THE  BOTTOM  OF  THE 
BAND,  THE  DENSITY  OF  STATES  CURVE  MUST  COINCIDE  WITH  THAT 
FOR  A  FREE  ELECTRON  BAND,  SINCE  ONE  OF  THE  CONDITIONS  ON 
EQUATION  (8)  WAS  THAT  II T  YIELD  THE  CORRECT  EFFECTIVE  MASS. 
AS  THE  ENERGY  INCREASES,  THE  PHYSICAL  CURVE  RISES  MORE 
SLOWLY  THAN  THE  FREE  ELECTRON  CURVE  DUE  TO  THE  EFFECT  OF 
THE  POSITIVE  TERMS.  BUT  Af  THE  ENERGY  RISES  STILL 
FURTHER,  THE  BAND  FLATTENS,  AND  THE  DENSITY  OF  STATES 
CURVE  CROSSES  THE  FREE  ELEOTROM  CURVE.  THIS  CROSS-OVER 
HAS  PROBABLY  OCCURRED  BEFORE  THE  FERMI  SURFACE  IS  REAOHEO, 
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SO  THAT  TH-  CFNSSTY  OF  STATCP  <N  THE  FlRMH  SURFACE 
SHOULD  3£  HIGHER  THAN  FOR  A  FREE  ELECTRON  GAS*  ONE 
WOULD  ALSO  EXPECT  DEVIATIONS  FROM  SPHERICA*  SYMMETRY. 

Beyond  the  Fermi  level,  the  density  of  states  should 

RISE  RAPICLT  TO  A  PEAK  AT  AN  ENERGY  !N  THE  NEIGHBORHOOD  OF 

the  Nj  level.  The  density  cf  statcs  will  then  fall*  At 

an  ENERGY  0.1|6  E.V.  AeO/E  THE  N  LEVEL,  THE  SECOND 

8ri;LLouin  Zone  begins  with  the  P-like  lfvel  N  *.  Thin 

6 

US  STILL  3/+  E.V.  BELOW  THE  TOP  OF  THE  FIRST  ZONE  AT 
H^,  80  THAT  THERE  IS  DC  NS  I  PER ABLE  OVERLAPPING  OF  THE 

bands.  The  total  bano  nioti-  from  r  to  is  I4..36  e„v., 

APPROXIMATELY  TWO  VOLTS  L'SS  THAN  WOULD  BE  EXPECTED  ON 
THE  BASIS  OF  ■*  HE  EFFEcriiVF  MASS  APPROXIMATION  AL<NE, 

The  occupied  i-orticn  of  the  band  has  a  width  of  about 
3.1  voltv  using  Ep  ,  E.j,  and  E^«  Equation  (8)  would  not 

li 

GIVE  A  SIGNIFICANTLY  DIFFERENT  RESULT. 

The  structure  of  the  x-ray  y  absorption  edge  in 

„  (28) 

POTASSIUM  HAS  BEEN  EXAMINED  BY  PlATT  .  DEVIATIONS 


28.  J »  Bo  Platt  ,  Phys.  Rev.  69.  377  (191,6). 

FROM  FREE  ELECTRON  ABSORPTION  IN  THE  DIRECTION  OF  LESS 
ABSORPTION  ARE  NOTICED  IN  A  REGION  BEGINNING  A  little 
MORE  THAN  1  E.V.  ABOVE  THE  FERMI  LEVEL.  ONE  WOULO 
EXPECT  0  STATES  TO  BE  IMPORTANT  IN  THIS  REGION,  AND  THESE 
MIGHT  SERVE  TO  REOUCE  THE  TRANSITION  PROBABILITY  TO  THE 


1 s  case  level. 
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IT  IS  INTERESTING  TO  COMPARE  THE  ENERQT  BANOS  OF 
POTASSIUM  WITH  THOSE  OF  THE  OTHER  ALKALI  METALS.  SlNOC 
RELIABLE  BAND  CALCULATIONS  HAVE  NOT  YET  BEEN  PERFORMED  FOR 
RUB  101 UM  ANO  OESIUM,  THE  PRINCIPAL  COMPARISON  WILL  BE  WITH 
sooium.  The  band  structure  of  sodium  has  been  stuoieo 
BY  MANY  WORKERS,  THE  MOST  RECENT  CALCULATION  BEING  THOSE 
OF  HAM^1^,  ANO  HOWARTH  ANO  JONES^^.  BOTH  CALCULATIONS 


29*  0.  J.  Howarth  ano  H.  Jones,  Proc.  Phvs. 

Soc.  65A,  355  0952). 

UTILIZE  THE  CELLULAR  METHOO,  BUT  WaM»S  WORK  IS  BASE0  ON 
the  Quantum  Defect  Method  ano  does  not  utilize  an  explicit 
POTENTIAL. 

Acoording  to  Howarth  ano  Jones  the  lowest  level 
at  the  Zone  corner  H  is  whioh  is  a  p  like  state. 

Next  oomes  and  then  H  •  (This  ordering  was  also 
FOUND  BY  Ham.)  H12  appears  to  be  the  lowest  state  at 
THIS  POINT  IN  POTASSIUM.  At  P,  THE  L0WE8T  STATE  IN 
800IUM  IS  PQ,  WHIOH  IS  8  LIKE.  IT  LIES  1*76  FT*  BELOW 
P;^  WHIOH  MIXES  P  ANO  0  STATES.  IN  POTASSIUM,  THIS  ORDER 
IS  RE VER8E0.  AT  N  THE  LOWEST  LEVEL  AOOOROINS  TO  HOWARTH 
ano  Jones  is  Ng,  but  Ham  has  founo  N(',  whioh  oontains  P 
FUNCTIONS  TO  LOE  LOWER  THAN  N  •  N(  APPEARS  TO  BE  LOWER 
IN  P0TAS8 I UM. 

It  is  likely  that  the  prinoipal  oifferenoes  in  band 

FORMS  BETWEEN  800IUM  ANO  POTASSIUM  CAN  BE  EXPLAINED  IN 
TERMS  OF  I NOREA8ED  IMPORTANCE  D  BANOS  IN  THE  LATTER  ELEMENT. 
IN  POTASSIUM,  THE  LOWE8T  STATES  WHIOH  OONTAIN  D  FUNOTIONS 


A. 


IT 
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TEND  CONS  ISTEN'LY  TO  LIE  LOWEP.  THAN  STATES  CONTAINING  ONLY  S 

OR  P  FUNCTIONS  (EXCEPT  OF  COURSE,  FOR  THE  BASIC  If.  LEVEL  T  ). 

8  I 

The  INCREA8E0  IMPORTANCE  OF  D  LIKE  LEVELS  SEEMS  ALSO  TO  LEAD 
TO  GREATER  DEVIATIONS  FROM  FREE  ELECTRON  LIKE  BEHAVIOR  IN 

this  case*  That  o  states  should  be  important  in  potassium 

WOULD  NATURALL'*  BE  EXPECTED  FROM  OBSERVATION  OF  THE  PERIODIC 
TABLE  OF  THE  ELEMENTS  SINCE  AS  SOON  AS  THE  ij  S  BAND  HAS 
FILLEO  IN  CALCIUM,  THE  0  BAND  BEG'NS  TO  BE  OCCUPIED  IN 
SCAN  0 1 UM »  l-OR  THIS  REASON,  CUE  WOULD  ALSO  EXPECT  TERMS 
IN  THE  E(k)  EXPANSION  TO  BE  MORE  IMPORTANT  IN  POTASSIUM 
THAN  IN  LITHIUM  OR  SOOIUM* 

VI.  HIGHER  BANDS 

IT  IS  INTERESTING  TC  OISCUSS  NOT  ONLY  THE  LOWEST  BAND 
BUT  ALSO  HIGHER  LEVELS,  IN  PARTICULAR  P  AND  D  LEVELS  TO  SEE 
TO  WHAT  EXTENT  THE  BANOS  ARE  SIMILAR  TO  THOSE  OF  OTHER 
ELEMENTS  WHICH  HAVE  THE  BOOT  CENTERED  CUBIC  LATTICE*  WE 
WILL  BE  OONCERNEO  HERE  P5IMARILT  WITH  THE  ORDER  OF  LEVELS 
WITHIN  A  GIVEN  8AND,  RATHER  THAN  RELATIVE  OROER  OF  THE  BANOS. 
WlGNER  HAS  GIVEN  AN  INTERESTING  ANALYSIS  OF  THE  STRUCTURE  OF 
P  AND  O  GAN  C  S  IN  THE  BOOY  CENTEREO  STRUCTURE  (JO),  BASED 


JO.  E.  WlGNER,  PROCEEDINGS  OF  THE  INTERNATIONAL 
Conference  of  Theoretical  Phtsic3,  Kisoto 
and  Tokyo,  September,  193J,  p«  650. 

ON  POSSIBILITIES  OP  SATISFYING  THE  BOUNDARY  CONDITIONS  WITHIN 
THE  CELLULAR  METHOO*  H’S  WORK  IS  IN  GENERAL  AGREEMENT  WITH 


THIS  CALCULATION 
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WlGNER  PREDICTS  FOR  THE 


P  BAND  THAT  THE  K®0  STATE 
l 


r 


15 


SHOULO  BE  HIGHEST,  THAT  THE  STATE  N(  SHOULD  BE  LOWEST 
WITH  A  LARGE  SEPARATION  OF  STATES  AT  THIS  POINT,  AND  THAT 

THE  STATE  H„c  SHOULO  BE  REASONABLY  CLOSE  TO  THE  BOTTOM 

15 

OF  THE  BAND.  THE  CALCULATION  CONFIRMS  THESE  PREDICTIONS. 

However,  Wisher  also  preoicts  that  should  lie  near  the 

MIDDLE  OF  THE  BAND  WHERE  HERE  IT  HAS  BEEN  FOUND  TO  BE  CON" 
S  I OERABLY  LOWER  THAN  HOWEVER,  SOME  OF  THIS  LOWERING 

MAY  BE  OUE  TO  THE  0  LIKE  FUNCTIONS  PRESENT  IN  Pj^. 


For  the  0  band  we  fino  N0  below  (V-  as  predicted. 

\  2  25 

States  of  it  perpend i cular  to  the  iii  directions;  i.e., 

AT  N  ARE  PREDICTED  TO  LIE  LOWEST.  SINCE  N(  MAY  BE  CON8IOEREO 
TO  BELONG  TO  THE  D  B AN  0  AS  WELL  A8  TO  THE  S  BANO,  THI8 
STATEMENT  IS  CORRECT.  NOTE  THAT  H^  ls  ALSO  LOW.  WiaNER'8 
0  I A6RAMS  OF  THE  RELATIONS  BETWEEN  P  AND  D  BANOS  ARE  IN" 
CORRECT  IN  THIS  CASE.  ALSO  THE  STATE  P^  LIES  BELOW  Ng. 

The  GENERAL  SHAPE  OF  THE  D  BAND  FOUND  HERE  IS  IN  AGREE- 

(17) 

MENT  WITH  A  RECENT  CALCULATION  FOR  BOOT  CENTERED  IRON'  1 ' 

BY  THE  O.P.W.  METHOO,  AND  IN  ALSO  WITH  THAT  FOUND  FOR  IRON 

(  31  ) 

BY  A  MOOIFIEO  TIGHT  BINOING  METHOD w  .  WE  OBSERVE  THAT 


31.  F.  Stern,  Ph.  0.  Thesis,  Princeton  University, 

1955  (unpublished). 

HERE  ALSO,  THE  SEPARATION  OF  THE  TRIPLY  AND  DOUBLY  DEGENERATE 
0  STATES  IS  MUCH  LESS  AT  THE  CENTER  OF  THE  ZONE  THAN  AT  THE 


2h 


CORNER  H,  ANO  THAT  THE  MAXI1MUM  SEPARATION  OF  THE  0  BAND 
IS  AT  THE  POINT  N»  EXCEPT  FOR  THE  FACT  THAT  THE  O.P.fl. 

IRON  CALCULATION  HAS  N|  ABOVE  N  ,  THE  RELATIVE  ORDER  OF 

4-  5 

THE  D  LEVELS  AT  A  3IVEN  SYMMETRY  POINT  IS  MUCH  THE  SAME 
IN  IRON  AMO  POTASSIUM.  WE  NOTE  THAT  IN  ALL  THREE  CASES  , 

THE  TRIPLY  DEGENERATE  D  BANO  STATE  IS  LOWEST  AT  K=0, 

BUT  THE  DOUBLY  OEGENERATE  STATE  IS  LOWEST  AT  H.  THE  ORDER 
OF  0  LEVELS  IN  POTASSIUM  IS  ALSO  IN  AGREEMENT  WITH  A  Ci‘ 

(52) 

LATION  FOR  BODY  CENTERED  TtTANIUM  ,  EXCi'T  IN  THAT  CASE, 

3 2.  B.  Schiff,  r^oc.  Phys.  Soc.  68a,  686  (1955)* 

A  MORE  EXTENSIVE  CALCULATION  WOULD  BE  DESIRABLE  IN  ORDER 
TO  MAKE  A  MORE  DETAILED  COMPARISON. 

There  is  a  simple  oevice  which  seems  to  reproduce  the 

CALCULATED  LEVEL  ORDER  WITHIN  A  GIVEN  BAND  VERY  WELL  IN 

this  case.  Let  represent  the  wave  function  for  the 

I  TH  IRREDUCIBLE  REPRESENTATION  OF  the  WAVE  VECTOR  K,. 

18  EXPANOEO  IN  0 R THO GON A L I  1 E 0  PLANE  WAVES. 

♦;  “  (9) 

WHERE  X„  .  „  18  GIVEN  0V  (3)«  THE  VECTOR  H  RUNS  OVER 

ALL  RECIPROCAL  LATTIOE  VEOTORS.  THE  LEVELS  MAY  BE  ORDERED 

t  2.  i 

ACCORDING  TO  THE  QUANTITY  (lK  *  H,  )  WHERE  (jj.  ♦  )  IS  THE 

O.P.W.  OF  LOWEST  ENERGY  BELONGING  TO  tK ' •  SlNOE  THIS  IS  AN 
ORDERING  AOOOROING  TO  KINE1IC  ENERGY,  THERE  IS  A  ROUGH 

(31  ) 

RESEMBLANCE  TO  THE  PRINCIPLE  OF  MAXIMUM  SMOOTHNESS w  . 


This  procedure;  says  in  essence,  that  the  order  or 

LEVELS  IN  A  GIVEN  9AND  IS  WHAT  lYOULO  BE  EXPEOTEO  FOR 
FREE  ELECTRON  BANDS.  IT  IS  NOT  SURPRISING  THAT  IT  GIVES 
THE  CORRECT  RESULTS  FOR  POTASSIUM.  IT  18  MORE  INTERESTING 
THAT  IT  SEEMS  TO  WORK  REASONABLY  WELL  IN  IRON.  EVEN  IN 
POTASSIUM,  HOWEVER,  ONE  MUST  NOT  SUPPOSE  THAT  THE  NUMERICAL 
ENERGY  VALUtS  OF  THE  LEVELS  WILL  BE  IN  GOOD  AGREEMENT  IV  I  TH 
A  FREE  ELECTS  *  Mr  rj,,  INSPECTION  OF  FIGURE  2  WILL  REVEAL 
THAT  SOME  LEVELS  W  E IL  CONNECT  El  V  BANOS  WHICH  MUST  DEPART 
SEVERELY  FROM  THE  FRIE  ELECTRON  FORM. 

The  SUCCESSFUL  COfPARtSON  OF  THESE  RESULTS  WITH  OTHER 
CALCULATIONS  FOR  THE  VARIOUS  BANDS  SUGGESTS  THAT  AT  LEAST 

THE  RELATIVE  OROER  OF  LEVELS  WITHIN  A  GIVEN  BAND  IS  NO! 

VERT  SENSITIVE  TO  DETAILS  CF  THE  POTENTIAL.  IT  IS  POSSI3LE, 
THEN,  THAT  THE  GENERAL  FORM!?  OF  INDIVIDUAL  S,  P  AND  0 
BANDS  ART  NOW  WELL  UNDERSTOOD  FOR  THE  BODY  CENTEREO  CUBIC 
LATTICE.  IT  18  IMPORTANT  TO  NOTE,  HOWEVER,  THAT  THE  SITUA¬ 
TION  IN  REGARD  TO  OVERLAP  CF  THESE  BANOS  IS  MUCH  LESS  CLEAR 

(  17) 

BECAUSE  THE  OVERLAP  IS  MUCH  MORE  SENSITIVE  TO  THE  POTENTIAL 
It  IS  INTERESTING  TO  NOTE  THAT  THERE  IS  APPARENTLY  A  VERY 
CC.N8  I  OERABLE  DEGREE  OF  OVERLAP  IN  POTASSIUM. 

Since  there  is  evioence  that  the  general  form  of  a 

BAAO  18  REASONABLY  INDEPENDENT  OF  POTENTIAL  IT  IS  INTERESTING 
TO  CON8IOER  THE  FORM  OF  THE  F  BAND.  AT  K*0,  THE  STATES  T 
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Tg8*  AMO  CONTAIN  F  FUNCTIONS.  AT  P  WE  HAVE  P^,  P^f 

AMO  P#>  AMO  AT  N  WE  HAVE  ^ * ,  Nf8,  Nj*.  THE  LOWEST 

PREDOMINANTLY  F  LIKE  STATE  I S  f ot ,  AMD  THE  TOP  OF  THE  BANC 

O' 

ALSO  APPEARS  TO  OCCUR  AT  T  WITH  f"2 * •  THIS  ORDER  IS  RE¬ 
VERSED  WITHS  H2’  LIES  BELOW  H gj..  We  f,no  ^  INTERMEDIATE 
BETWEEN  Toc  AND  I*  *  Ano  H. INTERMEDIATE  BETWEEN  H„  *  AND 
H^.  At  P#  THE  OROER  IS  APPARENTLY  P^,  P^,  ANO  Pj0  AT  N# 
N(‘  IS  APP''STNTLY  THE  LOWEST  FOLLOWED  BY  Ng 1 ,  Nj‘, 

N9*f  N^*#  **  O  Nj*  IN  THAT  ORDER.  WE  OON8IOER  HERE  THE 
8EC< NO  AND  THIRD  LEVELS  OF  THIS  CLASS  TO  BE  F  LIKE.  (NOT 
ALL  OF  THE  ENERGIES  OF  THE  LATTER  STATES  ARE  GIVEN  IN 

Table  III.) 


I  AM  INOEBTEO  TO  MR.  J.  SEGAL,  Mr.  E.  L.  HAASE,  ANO 
Mr.  D.  E.  Jack  for  assistance  with  the  rather  laborious 

COMPUTATIONS.  THE  ASSISTANCE  OF  PROFESSOR  R.  D.  WOODS  AND 

Mr.  S.  Berman  in  the  early  stages  of  this  calculation  is 

GRATEFULLY  ACKNOWLEDGED. 


TABLE  1 


The  quantity  rV  fs 

GIVEN  IN  ATOMIC  UN  1 rS  FOR  S,  P, 

AND  0 

ST* TES, 

AND  FOR  THE 

COULOMB  POTENTIAL. 

* 

rVs 

RVp 

RVo 

RVc 

.00 

58.00 

38,00 

58.00 

38.00 

.01 

36.73 

36.76 

36.70 

56.52 

.02 

35.61 

35.57 

35.49 

35.12 

.03 

34.51 

34.39 

34.30 

33.82 

*0  + 

33.46 

33.30 

33.19 

52.65 

.05 

52.53 

32.27 

32.15 

31-54 

gOS 

31.64 

31.52 

31.19 

50.54 

.07 

50.83 

30.41 

30.29 

29.61 

.03 

30.12 

29.58 

29.46 

28.75 

.09 

29.60 

28.77 

28.67 

27.93 

.10 

28.42 

28.02 

27.95 

27.15 

•  15 

24.78 

24.77 

24.70 

25.74 

.20 

22.22 

22.12 

21.92 

20.86 

.25 

20.10 

20.00 

19.62 

18.50 

.30 

13.53 

18.29 

17.67 

16.57 

•  53 

16.99 

16.82 

160  00 

15.02 

Jp 

16.08 

15.64 

14.64 

13.76 

.50 

13.41 

13.73 

12.56 

11.81 

.60 

10.91 

12.08 

11.06 

10.29 

.?o 

9.86 

10.48 

9.87 

8.98 

.03 

8.85 

8.84 

8.88 

7.83 

•  90 

7.93 

7.94 

8.02 

6.81 

1.00 

7.15 

7.17 

7.26 

5.95 

1.10 

6.47 

6.50 

6.56 

5.17 

1.20 

5.93 

5.95 

5.94 

4.54 

1  .Ip 

5.321 

5.093 

4.887 

3.598 

1.60 

4.262 

4.526 

4.054 

2.982 

1.30 

3.090 

3.462 

3.428 

2.595 

2. CO 

2.545 

2.839 

2.971 

2.357 

2.50 

2.286 

2.199 

2.550 

2.099 

3.00 

2.130 

2.049 

2.l4l 

2.025 

5-50 

2.055 

2.020 

2.076 

2.006 

J4..OO 

2.025 

2.009 

2.035 

2.001 

4.30 

2.009 

2.005 

2.013 

2.000 

t 


jr 
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TABLE  I  I 

Fourier  Coefficients  of  Potential  are  given  in  Rydbergs 

2  k  2 

as  Functions  of  the  Number  m  =  ( ■  -g— -  «  )  .  Only 

2 

Even  'ntegral  Values  of  m  Occur  for  the 
Sody  Centereo  Cubic  Structure. 


2 

M 

V8(m) 

VM) 

V0(M) 

Vc<M> 

0 

.8677 

.8692 

.871+8 

.7921+ 

2 

.28)4.1 

.2900 

.281+3 

.251+6 

k 

•231U 

.2369 

.2281 

.1965 

6 

.1865 

.1913 

.1821+ 

.1612 

8 

.1I+89 

.1525 

.11+39 

.1298 

10 

.121+9 

.1266 

.1192 

.1106 

12 

.111+3 

.111+5 

.1097 

.1035 

H+ 

.1000 

.0996 

.0973 

.0922 

16 

.0879 

.0675 

+0861+ 

.0821 

18 

.0775 

.0772 

.0769 

.0735 

20 

.  0686 

.0685 

.0688 

.0657 

22 

0651 

.0632 

.0636 

.0610 

2k 

.0585 

.0596 

.0591+ 

.0572 

26 

.051+0 

.0552 

.0562 

.0558 

TABLE  I  I  I 


Representation 


;iJEnav  States  in  Potajsi 

um  (Energies  in  Rydbergs 

of 

Determinant  No. 

OF 

Waves  Energies 

S-LIKE  STATES 

8 

1  55 

‘-e4304  +.8J46 

1). 

t>2 

+,2651  +.9780 

4 

ho 

+.0544  ♦‘3*  *Ul  5 

5 

P.2 

-.2239  +.2682 

P=!I.IKE  STATES 

4 

42 

+.6986  1.3277 

4 

:>3 

+.0621  +.9057 

U 

'Q 

-01262  +.9340 

4 

18 

-.1902  +.5561 

4 

20 

-.*.•769  .7813 

4 

20 

+*3561  1.0721 

O-LIKE  STATES 

4 

’  L  .8 

+.2200  +1.4814 

4 

54 

+•2640  +.9263 

4 

56 

+.5572  +1.2726 

4 

54 

-.1322  +1.1224 

3 

16 

+.4856  +1.3603 

4 

24 

+.0552  +.844-7 

4 

28 

+.6939  +1.2785 

4 

24 

+.4363  .8500 

F-LIKE  STATES 

3 

48 

+.2897  +1.8382 

2 

32 

+1.7401  4.1990 

3 

48 

1.4246  1.5285 

3 

56 

+.7046  +2.8097 

5 

36 

♦.5967  +1.3760 

3 

2'+ 

+.8945  2.0873 

TABLE  tV 


Lowes t  Gjge lvalue  by  Oroer  for  r  ( i n  Rydbergs).  The 

T  L3WEST  t  I  x  C3SHTH  ORDER  E  i1  G  K  N  V  A  l  U  E  S  ARE  ALSO 

Given. 

Lowest 

Oroer  £[\ 

1  -0.2371 

2  -.3834 

3  -.3900 

4  ~.kok9 

5  -.U125 

6  -.4232 

7  -.424.6 

a  -.4  304. 


Lowest  6  r.  1  qenvalues 

IN  EIGHTH  ORDER 

-.4304 
+.8346 
1 . 1478 
1.2621 
1.5309 
1.4525 


ft 


.£  ■'  \-.rrx 


CAPTIONS  FOR  -r I tittRES 


FIGURE  1, 

Brillouin  Zone  fo':  the  Body  Centered  Cubic  Lattice. 

FIGURE  2. 

Order  of  the  Energy  Levels  at  the  Four  Symmetry  Points 
r,  H,  P,  and  N.  All  Levels  with  E<1  My.  are  shown. 

FIGURE  3* 

Lowest  Energy  Band  along  100  Axis  in  the  Brillouin 
Zone  as  given  by  Equatiion  (8)  (Solio  Line),  and  as  given 
by  a  Free  Electron  Approximation  for  /  '  =  1.168  (Broken  Line). 

FIGURE  )+. 

Lowest  Energy  Band  along  110  Axis  in  the  Brillouin 
Zone  as  given  3v  Equation  (8)  (S^lio  Line),  and  as  given  sy 
a  Free  Electron  Approx «mat ton  for  s  1.168  (Broken  Line). 

FIGURE  5. 

Density  of  States  in  Arbitrary  Units  as  a  Function  of 


Energy 
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